Optimal scheme for estimating a pure qubit state via local measurements 
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We present the optimal scheme for estimating a pure qubit state by means of local measurements 
on TV identical copies. We give explicit examples for low N . For large N, we show that the fidelity 
saturates the collective measurement bound up to order 1/N. When the signal state lays on a 
meridian of the Bloch sphere, we show that this can be achieved without classical communication. 
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State estimation is a central topic in Quantum Me- 
chanics. Quantum measurements give only partial infor- 
mation about the state of the system under considera- 
tion. If an unlimited number of copies of a given state 
were available, one could in principle determine exactly 
this state provided an infinite number of measurements 
were performed. In the real world we have access to a 
limited number of copies and time for a finite number 
of measurements, and the best we can achieve is an es- 
timate of the state. It is, thus, of great importance to 
design optimal strategies which maximize the knowledge 
one can acquire about a general quantum state. 

Over the last few years a great deal of work has been 
devoted to the optimal estimation of states and many 
important features have been recognized [jl], [|, ||, [|, |[ |) , 
one of the most interesting ones being that collective 
measurements (CM) are more informative than individ- 
ual ones ^] . The experimental implementation of such 
measurements seems, though, quite involved. In practice, 
individual von Neumann measurements are far easier to 
perform. Moreover, one can show that optimal individ- 
ual measurements are of this type in the situations con- 
sidered in this letter [7|. The most general individual 
measurement procedures come under the name of LOCC 
(local operations and classical communication) schemes. 
In this framework, one allows a wide class of local oper- 
ations for which, depending on the outcome of the local 
measurement performed on a copy, appropriate transfor- 
mations can be applied on the subsequent copies of the 
state before measuring again. We consider the general 
situation in which the copies of the system may have 
never interacted in the past and may be macroscopically 
separated. Thus, classical communication becomes es- 
sential in these schemes. This is in contradistinction to 
collective measurements, for which the role of classical 
communication is played by the quantum entanglement 
of the measuring devices. In this LOCC context, some 
numerical simulations and experimental tests have been 
performed recently || [| . 

In this letter we will deal with pure states of a single 
qubit \ip). To each one of these states it corresponds a 
unique unit vector n on the Bloch sphere; i.e., = 
(1 + n ■ <t)/2, where a are the usual Pauli matrices. We 



will focus on two situations to which we will simply refer 
as 3D and 2D. In the former, no prior knowledge of n is 
assumed, whereas for the later (2D) n is known to lay 
on a meridian of the Bloch sphere. Our aim is to design 
LOCC measurements such that we can obtain the best 
estimate of ft and, hence, of We will investigate 

how good these LOCC measurements are as compared 
to the collective ones. For the latter, the mean fidelity F 
is commonly used as a figure of merit and many results 
can be found in the literature — in particular, the large 
iV asymptotic behavior of F is now known for the most 
interesting approaches |2|, ^] . 

For N identical copies of optimal CM lead to a 
fidelity that behaves as Fcm ~ 1 — 1/N for large N 
in 3D From JIo|, |ll| one can likewise compute the 
asymptotic fidelity in 2D, which is Fcm ~ 1 — l/(4iV). 
These results are the absolute upper bound for any mea- 
surement scheme, therefore within the LOCC framework 
the fidelity F can not approach unity at a rate larger than 
1 - 1/(4JV) in 2D (1-1/N in 3D) j|. In this letter we 
demonstrate that, rather unexpectedly, this asymptotic 
behavior can be achieved in 2D with just local measure- 
ments and no classical communication. In 3D classical 
communication seems necessary to saturate the asymp- 
totic CM bound, however, we have verified that for the 
optimal approach the fidelity reaches the CM regime very 
fast (for N ^ 12). Therefore, CM do not provide a sig- 
nificant improvement over local measurements. 

The estimation procedure goes as follows. After the N 
measurements (one on each copy) have been performed, 
a list of outcomes is obtained, which we symbolically de- 
note by x. Based on x, an estimate for \ip) can be guessed, 
\ipguess(x)) ■ The fidelity is then defined as the overlap 



| (V#guess(aO) I 2 = 



1 + n ■ M(x) 



fn(x), (1) 



where, as above, n is the unit vector on the Bloch sphere 
corresponding to the state and M(x) is that corre- 
sponding to IV'guess(^))- The average fidelity can be writ- 
ten as F = (/) = J2xl dnf n (x) P n (x), where P n {x) is 
the probability of getting the outcome x if the Bloch vec- 
tor is 71, and dn is the measure on the sphere in 3D (on 
the unit circle in 2D). Note that in both situations we 
have assumed no further prior knowledge of \tp). 



2 



Any local von Neumann measurement is represented 
by two projectors 0(±m) = (1 ± m • c?)/2, where m is a 
unit Bloch vector characterizing the measurement (in a 
spin system, e.g., m is the orientation of a Stern-Gerlach). 
We further note that for any unitary transformation U 
and for any m, there is always a vector in! such that 
0(ih') — UO(ih)U^ . So, any local operation on an indi- 
vidual copy of \tp) may be viewed as a redefinition of the 
Bloch vector that characterizes the measurement which 
will be performed on that copy. We hereafter adopt this 
point of view. For von Neumann measurements, the set 
of outcomes x can be expressed as a iV-digit binary num- 
ber x — i TV— 1 • • -Hii, were i k (= 0, 1) indicates that 
upon measuring on the fc'th copy, this is projected on 
the 0{{— Y k rhk] projection space. The most general local 
measurement is realized when we allow rfik+i to depend 
also on the list of previous outcomes ikik-i ■ ■ • *2*i = Xk 
(hence, x — Xn). We thus write ih(xu) instead of ink- 
Note that rh(xk) satisfy 



(2) 



The fidelity of a general LOCC scheme can be written as 



x=oo— o 



N 



fe=i 



where the last product is the probability P n (x). 

The optimal LOCC scheme is the one that maxi- 
mizes (J3J) over the Bloch vectors rh(xk) and M{x). Using 
the Schwartz inequality, it is straightforward to see that 
the best guess M{x) must be proportional to the vector 



V(x) = / dnn\[ 



1 + n ■ rh(xk) 



(4) 



fc=i 



i.e., M(x) — V(x)/\V(x)\. In this case, the maximum 
fidelity reads 



(5) 



For a fixed set of measurements {m(xk)} (optimal or not) 
and a given set of outcomes x, the guess (^) provides the 
best estimate of the signal state. This simple and general 
result does not seem to be conveyed in the literature. We 
next show how (^) can be used to improve the quality of 
state estimation schemes based on local measurements. 

Consider N = 2N copies of the state \ip) whose vector 
n is known to be on the equator (xy-v\ane) of the Bloch 
sphere (2D case). Since the expectation value of a is 
(i/j\a\ip) = n, the central limit theorem strongly suggests 
to adopt the following scheme. Let e\, e 2 be the two unit 
vectors pointing along the x and y axes respectively. For 
each i (= 1,2) perform N measurements of the observable 
<ji = 0(e*j) — 0(—ei). Assume we have obtained Na^ 



times the outcome +1 (consequently, N(l — m) times the 
outcome —1). One is driven to propose M oc J^i a i^i + 
— a i){~ 6j) (the mean value of these outcomes) as 
the Bloch vector of \if>). More precisely, 



Mi (a) 



20H-1 



E,(2^-i) 2 



(6) 



By doing so, the limiting behavior M (ip\a\ip) = it 

is ensured by the central limit theorem. With the tech- 
niques described below one can show that the fidelity for 
the guess (p) has the asymptotic expression 



Ft 



CL 



1 



3 1 

8iV 



(7) 



where CL stands for central limit. 

According to (^), however, this guess cannot be opti- 
mal. Let us show that indeed this scheme can be im- 
proved using (0) instead of (||). In this particular situ- 
ation Eqs. I and H read F OG = (1 + E Q l^(")l)/ 2 (OG 
stands for optimal guess) and 



*«> - n U, 



dn it 



x n (i±-) (1^) , (8) 

where we have used the shorthand notation a — (011, otz). 
To obtain the large N limit of Fog we use the formula 

( * ) ^(1 - g )d-°)>f -> 

\ a ^J 1 ' ^2^(1 - <z) 



x exp 



with q = (1 + rij)/2, and approximate the sum over the 
outcomes a by an integral. Performing the change of 
variables r ? - = 2a, — 1, we get 



d 2 r 



V 



d^r 



dn Qe 



(10) 



with E = J2 

r2 



n) 2 /(i 



and Q = P n = 



1 0* 

(N/27r)rj^ =1 (l — n\ ) -1 / 2 n. As it will become clear be- 
low, the terms of order l/v^N neglected in these approx- 
imations give no contribution to the final answer. The 
integral (|lfj) will now be evaluated with the aid of the 
saddle point technique. We use spherical coordinates and 
write it = (cos 0, sin 0), r = r (cos 7, sin 7). For large N, 
the integrand peaks at r = 1, which suggests to write 
r = 1 + £. We denote by <j) m or, equivalently, by n m , the 
location of the minimum of F, which can be computed as 
a power series in £ of the form <j> m = 7 + Ai£ + A 2 £ 2 + . . .. 
We next expand F and Q as a power series in t = 4> — 4>m 
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and write Q = L + Alt + Rt 2 + 



The exponential 



in (10) is seen to have the form 



2-E = e - 



- Si 



-xO(t J ) 



(11) 



where 2jB is the second derivative of E with respect to 
t, and the last exponential factor can also be expanded 
in powers of t. Neglecting contributions which vanish 
exponentially in N, one has 



V = / dnQe~^ E = 



V2ttJW 



L 



R 

JJB 



, (12) 



where dn = d<p /(2tt) = dt/(2n). At this point all the 
coefficients can be expanded in powers of £, e.g., B = 
B + Bi£ + R = Ro + Rx£ + . . ., and so on. Note 
however that the expansion of A starts at order £ 2 , i.e., 
A = A 2 e+A^ + . . .. For Q = nP = n m P+{n-fi m )P, 
we readily note that (O) yields 



V 



dn Pe 2 



I 



(13) 



where J/N is given by (|12j) with Q = (n—n rn )P. We have 
written a factor 1 /3\f in the last term to make explicit that 
this contribution is of order N -1 . Hence 



V 



dnPe~^ h + 



I 



O 



1 



(14) 



and we finally have up to 0(1/ "N 2 ) 

d 2 r\V\ = 1 - / d 2 r I dnP{\ - n ■ n m )e~^ E . (15) 



The leading term (unity) comes from the first integral 
in (|l4|), which is straightforward to compute exchang- 
ing the order of integration. Note also that the inte- 
gral of the leading term on the right hand side of (||), 
J d 2 rP exp(— 'NE/2), is one, and so is the sum over a on 
the left hand side independently of N. This shows that 
the terms of order l/\/N neglected in ( |Io|) indeed cancel. 
The subleading term in (|l5|) can be computed using the 
general expression (|l^) where now Q = (1 — n ■ n m )P = 
Rt 2 + 0(t 3 ) (i.e., L = M = 0). Moreover, only the lead- 
ing terms in powers of £ have to be retained, namely, A 2l 
Bq and i?o ( m particular, we just need Ai). This is so be- 
cause effectively 0{t 2 ) = 0(£ 2 ) = 0{l/N). The integral 
in (15) is then 



2tt 



diRo f°° , f XA 2 a 

, / arrexp< (r — 1) 

o ^WbZJo I 2 



2tt 



2tt 



d-y Rq 



(16) 



where again terms vanishing exponentially as N — * oo 
have been neglected. The cocmcicnt Ai is computed to 



be Ai = - cot 27. With this, B = 2, i? = (N/2w) esc 27 
and A 2 = 2 esc 2 27. The subleading term in ( [l5| ) yields 
— 1/ (4]Sf) , and substituting back in (pi) we finally obtain 



F< 



OG 



1 



1 1 

47V 



(17) 



(recall that N = 2N). As announced, this fidelity is 
larger than Fcl- Moreover, it saturates the absolute up- 
per bound given by CM even though classical communi- 
cation has not been used. In this sense, the behavior of 
the most basic local scheme (without classical commu- 
nication) in 2D is qualitatively similar to that obtained 
from CM provided the optimal guess (|j) is used. Working 
along the same line, we have computed the asymptotic 
fidelities Fql and Fqg m 3D. One has 



F, 



CL 



1 



6 1 

5iV 



F( 



OG 



1 



13 1 
12 N 



(18) 



where again we note that Fog > -Fcl- Despite this im- 
provement, the CM bound Fqm = 1 — 1/^V" is not satu- 
rated, although the subleading term of Fog is only 8% 
less than that of Fcm ■ 

So far classical communication has not been exploited, 
i.e., the Bloch vectors rn{xk) were both x-independent 
and non-optimal. We now turn to the full-fleshed LOCC 
schemes. Hereafter only the general case 3D will be con- 
sidered. 

The first non trivial case is N — 2. Here, x takes four 
possible values, 00, 01, 10 and 11. There are three inde- 
pendent vectors, namely, m(0), m(00), m(01) (the other 
three are obtained using Eq. g). The first vector, m(0), 
is arbitrary and we take m(0) = z. The optimal fidelity 
is obtained by maximizing (J5j) with respect to m(00) and 
rn(01). A straightforward calculation gives 



fc=0,l 



£ !?(*)! = £ £ (|sin^| + |cos% 



#0fc I , I #0fc 



(19) 



where 9 X is the polar angle of the vector m(x). The 
maximal value of ( |l9|) is attained for 6? 00 = #01 = tt/2. 
Notice that 9 00 and 6 m are maximized independently, 
so m(OO), m(01) do not need to be equal; they are only 
required to be orthogonal to m(0). Substituting back 
in f| we find F^ = (3 + V2)/6 (see also (|). From 
Eq.^ we obtain the optimal guess 



M {2 \x) 



m(x 2 ) + m(xi) 
V2 



s(x), 



(20) 



hence, e.g., Af( 2 )(01) = [m(01) + m{l)]/V2 = [rn(01) - 
m(0)]/\/2. 

The case N = 3 is very similar. The optimal Bloch 
vectors, rh{x{), rh(x2), rn{x^), are found to be mutually 
orthogonal. Since there is no further constrain, one can 
choose the three fixed (i.e., independent of x) vectors 
w(xfc) = e*fc, k = 1,2,3. This shows that for N = 3 
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(as well as for N — 2) the optimal estimation schemes 
based on local measurements do not require classical 
communication. For each outcome x the optimal guess 
is M^(x) — {mfas) + m(x 2 ) + m(xi)]/V3, which is 
a straightforward generalization of M^- 2 \x), and yields 
— (3 + v / 3)/6. These results could somehow be an- 
ticipated: if 0(m)\ijj) ^ we can only be sure that the 
Bloch vector of is not —in. Intuition suggest to use 
the subsequent copies of \ip) to explore the plane orthog- 
onal to m. Thus, the optimal Bloch vectors rh(xk) tend 
to be mutually orthogonal. 

The case ./V = 4 is more complex, since four mutu- 
ally orthogonal vectors cannot fit onto the Bloch sphere. 
The solution exhibits some interesting features. First, 
the optimal Bloch vectors now depend on the outcomes 
of the previous measurements. Therefore classical com- 
munication does play a crucial role for N > 3. However, 
fh{xi) _L m(x2) and, as before, one can choose rh(xi) — 
e%, for i = 1,2. Only for the third and fourth measure- 
ment one really has to take different choices in accordance 
to the sequence of the preceding outcomes. The Bloch 
vectors of the third measurement can be parametrized 
by a single angle a as rh(x 3 ) = cosaui(x) + smavi(x), 
where U\(x) = m(xi) x m{xi) and V\(x) = U\{x) x s(x) 
and s(x) is defined in (|2(|). The optimal value of this 
angle is a = 0.502. We cannot give any insight as to why 
this value is optimal. However, in agreement with our in- 
tuition, m{x^) _L s{x), i.e., the third measurement probes 
the plane orthogonal to the Bloch vector one would guess 
from the first two outcomes. Two angles are required to 
parametrize the vectors of the fourth measurement. They 
are given by m[xi) = cos 71*2(2;) + sin 7 1/2 (a;), where 
U2(x) = s(x) x fh(xs), V2(x) = cos j3 rh{x^) — sm.(3s{x). 
The optimal values of these angles are (3 = 0.584, 7 = 
0.538, and the corresponding fidelity is — 0.8206. 
This is just 1.5% lower than the absolute bound Fq^ = 
5/6 = 0.8333 attained with CM. We also give the values 
of the maximal LOCC fidelities for N = 5,6. They are 
F( 5 ~> = 0.8450 and F< 6) = 0.8637. It is interesting to note 
that for N > 3, it pays to relax optimality at each step. 
Hence, one-step adaptive schemes || [| are not optimal, 
though the differences are very small; e.g., for N = 4, 

> ^idiptivc = (15 + V9l)/30 « 0.8180. 
Having learnt from the low N cases, we are in the 
right position to compute the asymptotic fidelity of 
this scheme. For that, we take inspiration in varia- 
tional methods as follows. Suppose we have performed 
a large number Ao = y/N of measurements and ob- 
tained the guess Mq. It is clear that the subsequent 
27V = N — No guesses will hardly differ from Mq. We 
hence substitute in (||) the ansatz M(x) w Mo cos to + 
sin uj(u cost + usinr), where u, v are two unit vectors 
which along with Mq form an orthogonal bases, u> — 
X^(2a u - l) 2 + (2a„ - l) 2 , tanr = (2a v - l)/(2a u - 1), 
and A is a variational parameter. As above, Na u (Na v ) 



is the number of times we obtain the outcome +1 when 
we measure a ■ u (<r • v) . Note that in average u) will be 
small since we expect a u ^ v w 1/2, and we need to retain 
terms up to order ui 2 . Putting all this together one gets 
from (||) 

F>l- ( l-A)*(l-Fo)-^ + ..., (21) 

where Fq is the optimal fidelity for Nq measurements and 
the dots stand for subleading terms in inverse powers of 
and ATq. We readily see that the optimal choice is 
A = 1, which leads to F m 1 - 1/JV [|||. Hence, our 
LOCC scheme does saturate the CM bound. Further- 
more, numerical analysis reveals that the CM regime is 
reached for values of N as low as 12. 

In summary, we have obtained the optimal LOCC 
estimating scheme for general qubit (pure) states and 
shown that its fidelity saturates the collective measure- 
ment bound. For states that are known to lay on a merid- 
ian of the Bloch sphere (2D case) we have explicitly given 
a scheme whose fidelity saturates this bound without in- 
voking classical communication. 
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